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Low energy excitations of the kagome antiferromagnet and the 
spin gap issue. 
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Abstract. - In this paper we report the latest results of exact diagonalizations of SU(2) invariant 
models on various lattices (square, triangular, hexagonal, checkerboard and kagome lattices) and 
revisit critically in this light some of our previous analysis. We focus on the low lying levels in each 
S sector. The differences in behavior between gapless systems and gapped ones are exhibited. The 
plausibility of a gapless spin liquid in the Heisenberg model on the kagome lattice is discussed. 
A rough estimate of the spin susceptibility in such an hypothesis is given. The evolution of the 
intra-5* channel spectra under the effect of a small perturbation is consistent with the proximity 
of a quantum critical point. We reemphasize that the very small intra- S channel energy scale 
observed in exact spectra is a very interesting information to understand the low T dynamics of 
this model. 
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The RVB concept of a quantum Spin Liquid due to 
P.W. Anderson in the beginning of the seventies has mo- 
tivated a very large number of works both theoretical and 
experimental [1,2]. Very recently a specifically interest- 
ing candidate for a spin liquid has emerged: the Herber- 
smithite characterized by a fluctuating behavior down to 
temperatures of the order of 10 -4 times the coupling con- 
stant [3-7]. In this compound the magnetic Cu ions are 
located at the vertices of a kagome lattice and the interac- 
tions between these spin-1/2 ions arc supposed to be fairly 
isotropic. The possibility that a simple Heisenberg model 
might capture the physics of this compound, has there- 
fore attracted interest, although a significant amount of 
defects and other interactions, e.g. Dzyaloshinskii-Moriya 
interactions, are likely to be present. However, in spite 
of many studies, the nature of the ground-state and first 
excitations of the spin-1/2 kagome Heisenberg antiferro- 
magnet (KHA) is still disputed [8-16]. 

Exact diagonalizations [8-10, 17, 18] have very early 
shown that something special happens in the KHA. n- 
points correlations functions (up to n=8) decrease at short 
distances [8,16-21]. The spectral properties observed on 
these small samples (up to N = 36 sites) are quite un- 
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usual: these spectra exhibit two low-energy scales: be- 
sides a finite-size spin-gap scale, there is a second much 
lower energy scale inside each sector of total spin S. A 
crude 1 /N extrapolation of the finite-size spin-gaps to the 
thermodynamic limit lead to an "estimate" of a putative 
thermodynamic spin-gap of about 1/20 of the coupling 
constant [10]: an analysis and a result that we want to re- 
visit critically at the light of recent more extensive studies 
of the spectra of samples up to N = 36 spins. 

Many other numerical approaches have been used to 
attack this problem: CORE renormalization approaches 
(Budnik et al. [11] and Capponi et al. [22]), sophisticated 
series expansion (Singh et al. [12,15]), DMRG calculations 
(Jiang et al. [16]), and more recently MERA renormaliza- 
tion procedure (Evenbly et al. [23]), give evaluations of 
the ground-state energy and for most of them of the spin 
gap. 

Since 1998 and our publication [10] the performances of 
computers have increased noticeably (not enough to com- 
pute the spectrum of samples much larger than 36) but 
sufficiently for computing a large number of levels in any 
S sector of various clusters up to JV = 36 spins and thus 
allowing a full analysis of the exact T = magnetism 
of this "large molecule". The object of this paper is to 
display these results and compare them to the behaviour 
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of the Heisenberg antiferromagnet on various other lat- 
tices: square (SHA), hexagonal (HHA), triangular (THA), 
checkerboard (CHA) of similar sizes. We show that up to 
this size, the KHA has a magnetic behavior which com- 
pares extremely well to other antiferromagnets known to 
be gapless in the thermodynamic limit and that, the ex- 
trapolation of such finite size data although it does not 
rule out a finite spin-gap, cannot yield a significant mea- 
surement of an eventual very small spin gap: our 98' paper 
was on this point misleading [10]. From the methodologi- 
cal point of view the analysis of the data that we perform 
in this paper can readily be applied to any results obtained 
by other numerical finite size calculations (DMRG etc.). 
The third message of this paper is that the much lower 
energy scale that appears inside each 5 sector play a key 
role in the true dynamics of this magnet and suggests a 
critical or quasi-critical behavior. 

Gapless and Gapped magnets through exact diagonaliza- 
tions: ground-state energy versus total spin 

In Fig. 1, we display the lowest eigen- levels of the SHA, 
THA, HHA, KHA and CHA samples of 36 spins (with 
periodic boundary conditions) and study how they evolve 
with the total spin value 5 of the sample. The hamiltonian 
is 

H = Y,3i- *i (!) 

(id) 

where s*j is the individual spin 1/2 at site i and the sum 
runs over pairs of nearest neighbors (5 = Sj). 

In the checkerboard lattice, also known as the py- 
rochlorc slab lattice, tetraedra share corners and are dis- 
posed on a checkerboard: all couplings along tetraedra 
edges are identical. 

Fig. 1 top shows that in most of the cases (except CHA), 
the evolution of the first level in each 5 sector is extremely 
well described by a 5(5 + 1) behavior which extends to 
5 = S m ax = N/2 (not shown in the figure for the sake of 
clarity) 1 . 

Such a behavior can be described by the equation: 



e N (S) = e N (0) + a N 



5(5 + 1) 

N 2 



(2) 



where ej\r(5) is the lowest energy per spin in the 5 channel 
for a sample of size N. The limit of Eq. 2, when the size of 
the sample goes to infinity is the well-known formula of the 
energy versus magnetization of classical antiferromagnets: 



e(m) = e(0) + — m 2 
2x 



(3) 



where m is the magnetization (0 < m < 1) of the sample 
and x the dimcnsionlcss susceptibility. (With this defini- 
tion, the magnetic moment per spin in a magnetic field H 
reads fj, = {g^ B /^) 2 xH and a N = 2/xn)- 



1 In fact there is a tiny cusp atS = 6(l/3of the maximum mag- 
netization) both for the triangular system and the kagome system 
which signals a magnetization plateau, and a change of behavior be- 
tween Smax — 1 and Smax which is the manifestation of localized 
states in the nearly magnetized KHA [24] . 




Fig. 1: (Color online) Lowest energy per spin of the nearest 
neighbor Heisenberg model on different lattices versus total 
spin: square (green squares), triangular (bluetriangles) hexag- 
onal (pink diamonds) 2-dimensional pyrochlore slab(open blue 
squares) and kagome (red bullets). Top: exact energies are 
displayed versus S(S + 1). The full lines are linear regressions 
according to Eq. 2. For the pyrochlore lattice (dot-dashed line) 
the linear regression only involves the eigenstates with total 
spins larger or equal to 4. Bottom: the eigen-levels are dis- 
played versus 5*. The pyrochlore results are linear in S for 5* 
ranging from to 3. All samples have 36 sites and the full sym- 
metry of the infinite lattice (periodic boundary conditions). 



The low energy effective behavior described by Eq. 2 is 
expected from first principles for a uniaxial ordered anti- 
ferromagnet [25], as is the case of the SHA and HHA on 
the square and hexagonal lattices [26] . In such a case, 
Eq. 2 describes the free precession of the order parameter 
of the magnet (which has the symmetry of a rigid rota- 
tor) in an SU(2) invariant environment. A slightly more 
complicated behavior is expected for a biaxial magnet like 
the THA (on the triangular lattice): the order parameter 
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Fig. 2: (Color online) Top: Evolution of the ground-state en- 
ergy per bond with the sample size (these valuesare extracted 
from the fit of the low lying levels of the different spectra to 
Eq. 2. The variation of the energy with size and shapes of sam- 
ples is rather smooth (in spite of the fact that many samples 
used for this graph have not all the symmetries of the infinite 
lattice). Bottom: evolution of the spin susceptibility with the 
sample size. The size effects are in first approximation propor- 
tional to 1/TV 3//2 for the ground-state energy per spin and to 
I/TV 1/2 for the spin susceptibility. 



has now the symmetry of a quantum top and two different 
susceptibilities are needed to describe its free motion [27]. 
It can indeed be observed in Fig. 1 that the fit of the low 
lying levels of the THA with Eq. 2 is not as good as in 
the case of collinear Neel order (SHA et HHA) . In fact the 
Pearson product-moment correlation coefficient of the fit 
of the 6 first levels (Stot = • • ■ 5) of the exact spectra to 
Eq. 2 is 1 for HHA, resp. 0.99997 for SHA, 0.99996 for 
KHA and only 0.99849 for THA. 

In addition to the above discussed data, Fig.l-top ex- 



hibits the low lying levels of the spectrum of the CHA 
(a two dimensional pyrochlorc lattice). This system is 
a good example of a Valence Bond Crystal with a large 
gap of the order of 0.7 [28,29]. The S dependance of the 
CHA energy is clearly different from that of the previous 
models. For high values of the total spin CHA energy dis- 
plays the standard semi-classical quadratic behavior de- 
scribed in Eq. 2. But for S < 4 the energy of the magnet 
deviates clearly downwards from this quadratic behavior, 
and appears dominated by a linear in S behavior (Fig.l- 
bottom) 2 . The following heuristic law is then a good first 
order approximation of the small S behavior: 3 



E N (S) = E N {0) + A N S. 



(4) 



where Ew(S) is the total energy of the ground-state of the 
TV-sites sample in the S channel. In a magnetic field this 
equation becomes 



E N (S) = E N {0) + A N S z - g f i B HS z 



(5) 



where Ajv, the spin gap (~ 0.7 in the case of the CHA) is 
a measure of the critical field H c needed to observe a mag- 
netic answer of the system (in units where g hb = 1)- At 
that point we can do the contact between the usual point 
of view in term of one-particle excitations and our many 
body approach. The so-called Bose-Einstein condensation 
of spin excitations at H c translates in our approach in a 
collapse of a finite fraction of cigen-lcvels with different to- 
tal S z when H = H c = A. In our 36-sitc CHA sample, the 
linear in S behavior of Eq. 4 implies that the eigen-levels 
with total spin components S z — 1, 2, 3 (which represent 
respectively states with 1, 2, 3 bosons) become degenerate 
when H reaches H c . 

The perfect adequation of Eq. 2 for the description of 
the low lying levels of the KHA (for all sizes up to TV = 36) 
opens a question: is it legitimate to ascribe the finite-size 
gap of the spectrum of this sample to a true spin gap as it 
can be perfectly described in an alternate, ultimately gap- 
less picture? Either a saturation of the gap with finite 
size 4 or a magnetization curve, as described in Eq. 4, 
should be observed to ascertain a true thermodynamic 
spin-gap. With the knowledge of the exact spectra of the 
KHA up to TV = 36, it seems impossible to decide if KHA 
is a gapless system or a gapped one with an extremely 
small gap. 



2 The downwards deviation is the manifestation that resonances 
which stabilize the quantum magnets arc larger and larger when S 
decreases, as expected. 

3 Indeed, specifically in one dimension, this mean-field approach 
and Eq. 4 is certainly a bit naive. Non analytic terms dominate the 
behavior of the spin susceptibility just above the spin gap, but the 
overall thermodynamic behavior is correct 

4 In fact neither the Contractor Renormalization approach [11] 
(CORE), nor the Density Matrix Renormalisation Group [16] 
(DMRG) do observe a real saturation in the gap with increasing 
sizes up to 81 (resp. 108). In that last approach this saturation is 
only obtained for quasi uni-dimensional samples with 3 or 4 legs. 
But the same saturation has not been exhibited for larger bands. 
The MERA renormalization [23] approach may be promising but 
has not given up to now an estimate of the gaps. 
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In the hypothesis of a gapless system with linear in k 
low lying excitations, e(0) (resp. x) would scale as N~ 3 ^ 2 
(resp. N" 1 / 2 ). The fits in Fig. 2 are not contradictory 
with this hypothesis. If this picture is valid up to the 
thermodynamic limit, the extrapolated value of the spin 
susceptibility is lim^ —,00 (xn) ~ 0.09 ± 0.015. Within 
the error bars (which arc rather large), the extrapolated 
value of the susceptibility compares favorably to the bulk 
susceptibility measurement of the Herbcrtsmithitc [6] and 
is consistent with the high temperature series data [30]. 

The large finite-size spin gap appears in this point of 
view as a simple manifestation of the finite value of the 
uniform spin susceptibility and of the spin quantization. 
In this aspect the behavior of KHA is not very different 
from that of SHA, HHA or THA. But, as we will show be- 
low, the situation is quite different for the intra-5-channcl 
excitations, which do not suffer of the same limitation due 
to quantization and allow the exploration of a much lower 
energy scale. 

Neel ordered magnets, energy scale of the excitations 
In the two first situations SHA and THA (Fig. 3 a and 
b), all the low lying levels in each 5 sector above the 
ground-state describe the spin waves. The analysis of the 
5=1 sector is specially simple 5 : one easily recognizes 
the soft mode(s) of the SHA (point X of the first Brillouin 
Zone (BZ) Fig. 3 a) and of the THA (points V and K of 
the BZ of the THA Fig. 3 b). For the THA the dispersion 
curve can directly be compared to the most sophisticated 
series calculation [31]. The agreement is qualitatively total 
(flat zone, local minima) and quantitatively rather good. 
In spite of the small size of the samples, the order of mag- 
nitude of the spin-wave bandwidth (2 for the SHA and 
only 1 for the THA) as well as the general behavior of the 
dispersion curve are in correct agreement with thermody- 
namic calculations. Note that in the typical bandwidth of 
one magnon excitations, the exact spectra in the 5=1 
sector exactly display one eigen-level per wave number in 
the SHA case and approximately the same number in the 
triangular case (crosses in Fig. 3 a and b). (In that latter 
case the bi-axial nature of the magnet brings some minor 
complications). 

Specificity of the kagome antiferromagnet 
By comparison the spectrum of excitations of the KHA 
Fig. 4 is qualitatively and quantitatively totally different 
from the two previous systems. It exhibits a consider- 
able density of low lying excitations, seemingly continuous 
above the ground-state in each S sector. At same sizes, 
the finite size gaps inside an 5 sector are much smaller. A 
molecule of 36 spins with Heisenbcrg nearest neighbor cou- 
plings and the coordination symmetry of a square would 
have internal Bohr frequencies for AS = transitions in 
the 5 = 1 sector of the order of 1.5. This number goes 



J The analysis in the 5 = sector is more involved, selection rules 
on the coupling of angular momenta imply that the excited levels in 
the S = sectors are uniquely two-magnons excitations when the 
order parameter is a rigid rotator, whereas for a top- as is the case 
of the THA- there are both one and two magnons excitations. 
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Fig. 3: (color online) Energy above the ground-state of the 
lowest energy levels, as a function of wave vectors, in the lowest 
spin sectors S for the N — 36 clusters of the square (top) and 
triangular lattices (bottom). Insets show the location of the 
wave vectors in the Brillouin zone. These clusters have all the 
symmetries of the infinite lattice. The full spin wave dispersion 
of the square and triangular antiferromagnet is indicated by the 
blue-red crosses; the soft modes by the red circles. 



down to 0.5 for the triangular coordination and to 3. 10 -3 
for the kagome network. This frequency scale at this size 
is only slightly larger than the one observed in Herbert- 
smithite -in which a fully developed dynamics down to 
temperatures of the order of 2.5 10~ 4 [4,7] has been ob- 
served. The similar low energy of the AS = transitions 
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Fig. 4: (color online) Same as 3 for the kagome lattice. Note 
the differences in the vertical scales for these three samples. 
The density of low lying states for the kagome sample is a few 
hundred time larger than in the two other magnets of the same 
size. The convergence of the Lanczos process is thus more 
tedious and is only achieved below the vertical red lines. In 
the range of non-converged levels, the displayed spectrum is 
less dense than in reality. 

in the S = sector equally indicate fluctuation of singlets 
down to very low temperatures. A measurement of the en- 
ergy scale inside each S sector on larger samples would be 
of the utmost interest as it might match the experimental 
scale observed in Hcrbertsmithite for N ~ 108. Presently 
impossible with ED it should be reachable in DMRG [16] 
or MERA [23]. With some caveat, analysis in the dimer 
basis could also be interesting [34,35]. 

Proximity to a quantum critical point? 

The huge density of low lying states in each S sector, 
and its evolution with the system size [10] which is the 
specificity of KHA, can be interpreted as the absence of 
an intrinsic energy scale (at least in these small samples). 
This opens the question of the criticality [13,14,32]. We 
show in Fig. 5 the evolution of the low lying levels in 
the S = (resp. S = 1/2) sector of the N = 36 (resp. 
N = 27) sample under the action of an anti-ferromagnetic 
second neighbor coupling. Under the action of this per- 
turbation, a gap opens linearly just above the ground- 
state. The levels that would embody a VBC symmetry 
breaking [12, 15, 33] arc immediately destabilized as well 
as all other excited states. The ground-state (resp. in 
the S = 1/2 and 5 = channels) has all the symmetries 
expected for the g.-s. of the q = Neel order. At this 
scale and inside each S sector KAH looks like as it was 
critical 6 . In refs. [9,37], we concluded that a non zero 

6 The same situation was observed when destroying the three sub- 



l il^ ll ± llll i llll ± llll l l 



0.4 - = s = 



0.3 - = _Z 



W - = s - 

CO 0.2 - " = " 

tafl - z - z 



0.1 - Z E - 

; H: - (a) ; 

: ~.= n=27: 

I-, T ~ 

0.05 0.1 0.15 0.2 









| I 


' § 1 










* 


0.2 


















e 










0.15 






- 

s 




0.1 




— 9? 

— * 

= 5 






0.05 


- e 
_ ^ 

- 8 

- e 


i g 




(b) " 

N = 36; 





A 






, x , 







0.05 


0.1 



Fig. 5: (color online) Evolution of the low lying levels of the 
spectrum of the Ji — J2 model on the kagome lattice versus J2 
(Ji = 1) for a symmetric N — 27 sample (a) and for N = 36 
(b). The total energy of the eigen- levels is measured rela- 
tive to the absolute ground-state for each value of J2 (black 
half circles). Note that a gap opens quasi linearly just above 
the ground-state. The levels that would embody a VBC-1 or- 
der [12, 15, 33] are distinguished by black circles (symmetric 
pinwheel) and red crosses (antisymmetric one). 

critical coupling (cither a second neighbor interaction [9], 
or Dzyaloshinski-Moriya [37]) was needed to achieve full 
q = Neel order. These two determinations are implicitly 
dependant of the large energy scale set by the finite size 
gap (Aat =36 = 0.2). At the light of the present study on 
can at best ascertain that the values obtained in these pa- 
pers are certainly upper-bounds. The present study can- 
not cither ascertain that the system is critical. In some 
sense it may be instructive to compare the finite-size sit- 



lattice long range order on a triangular lattice by a cyclic 4-spin 
exchange [36]. 
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uation to a non zero temperature regime: in the vicin- 
ity of a T = critical point, it is widely accepted that 
there exists a T / quantum critical regime which ex- 
tents on a finite range of couplings around the quantum 
critical point. It seems quite reasonable to argue that the 
behavior observed in rcfs [9, 37] around the pure Hcisen- 
bcrg kagome point looks like a quantum critical regime. 
Such an observation does not preclude at lower tempera- 
tures/larger sizes a subsequent crystallization of the pure 
KAH in a very large Valence Bond Crystal or a small gap 
RVB phase, but it gives some support to the idea that 
the system could be near a critical point, and that in a 
non negligible range of energy, time and length scale it 
could behave as " quasi-critical" . Such a point of view is 
also supported by the conclusion of a finite temperature 
DMFT calculation done some years ago [38] . 

In this paper we have explicitly shown that the KHA 
spin-gap measured in numerics for sizes up to 36, could be 
an extrinsic property: on such small samples it is impossi- 
ble to distinguish between a gapless system and a system 
with a very small gap. We have shown that the intra-5 
channel spectra could signal the proximity of a quantum 
critical point. We suggest that polarized inelastic neutron 
scattering on Herbertsmithite at temperature low enough 
but larger than the residual Dzyaloshinski-Moriya pertur- 
bation would bring an interesting insight on the spin-1/2 
compound, with a possible comparison to SrCrGaO [39]. 

* * * 
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